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Abstract. The paper presents a reformulation of some of the most basic 
entities and equations of linear elasticity — the stress and strain tensor, the 
Cauchy Navier equilibrium equations, material equations for linear isotropic 
bodies - in a modern differential geometric language using differential forms 
and lie derivatives. Similar steps have been done successfully in general rel- 
ativity, quantum physics and electrodynamics and are of great use in those 
fields. In Elasticity Theory, however, such a modern differential geometric ap- 
proach is much less common. Furthermore, existing reformulations demand a 
vast knowledge of differential geometry, including nonstandard entities such as 
vector valued differential forms and the like. This paper presents a less general 
but more easily accessible approach to using modern differential geometry in 
elasticity theory than those published up to now. 



1. Introduction 

Modern differential geometry often allows for a clearer, more geometric approach 
to solve physical problems than Gibbsian vector calculus or Ricci calculus (i.e. the 
index-based form) of Riemannian geometry. Even more importantly, equations 
stated in modern differential geometric terms are not bound to the use of special 
classes of coordinate systems, such as cartesian or orthonormal ones. It provides 
useful and powerful concepts such as differential forms and Lie derivatives and 
with them a generalization of vector calculus expressions like rotation, gradient 
and divergence and of various integral theorems like Gauss' divergence theorem 
or Stoke's integral theorem on planes to just one integral theorem. Furthermore, 
coordinate transformations in curvilinear coordi-nate systems are often simplified. 

In general, it is desirable to express equations in a connection and metric ten-sor 
free manner, using only differential forms and Lie derivatives, thus simplifying co- 
ordinate system changes. Unfortunately it is not always possible to express vector 
or tensor equations as equations in such a way. One reason of this is the anti- 
symmetric nature of differential forms of degree >2. While e.g. the gradient of a 
scalar function and the rotation and divergence of a vector field can be expressed 
using differential forms, this is not possible with the gradient of a vector field . It 
is, however, possible to reformulate certain differential equations of vector fields as 
ones of differential forms. 

In electrodynamics, differential forms are already a quite popular way to refor- 
mulate maxwell's equations etc. |Sel97| In continuum mechanics, however, there 
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are only few steps towards using differential forms, most notably from J. E. Mars- 
den |Mar94| and V.I. Arnold |Arn89| . In |Kan07| . the stress and strain tensors 
are introduced as vector or covector valued differential forms of various degree. In 
that approach, the physical meaning of those tensors are clearer than in the tradi- 
tional way of expressing them interchangeably as contravariant, covariant or mixed 
tensors. However, vector or covector valued differential forms are a nonstandard 
subject in differential geometry and usually do not simplify calculations. 

For that reason, in this paper a path between traditional tensor analysis and the 
very general and mathematically rigorous approach of fMar94j is taken, thereby 
restricting to linear elasticity. Due to their relevance, only static problems without 
body forces are considered. 

Absolutely necessary to determine deformations inside a body is a suitable sys- 
tem of differential equations for displacements, stresses or strains, respectively. 
Since displacement fields can be considered as fundamental in the sense that stress 
or strain fields can be determined from them, the cauchy-navier displacement equa- 
tions shall be reformulated in chapter ([2]). The resulting set of equations have only 
the displacement variables as unknowns, so displacement boundary conditions are 
easy to handle in this system of equations. In order to allow also traction boundary 
conditions and to relate stresses, strains and displacements, it is necessary to ex- 
press the well-known stress and strain tensor of linear elasticity in the same modern 
differential geometric formalism . This will be done in chapter ([3]). 

It shall be noted that throughout the paper it was assumed that manifolds have 
sufficiently smooth boundaries and boundary values are also sufficiently smooth. 

2. Cauchy Navier's Displacement Equations of Static Linear 

Elasticity 

2.1. Basic Definitions. Problems in static linear elasticity consist of finding ei- 
ther stress fields that satisfy the equilibrium equation V • ^ = or, equivalently, 
displacements that fulfil the Cauchy Navier equation |Tim70j inside a suitable sub- 
manifold B oi E"^: 



In this article, we will regard the latter equation as fundamental because stresses 
and strains can easily be calculated from displacements. Also, there is no need for 
additional compatibility equations to be satisfied as is the case for the equilibrium 
equation for stresses or strains. 

Displacements in elasticity theory are naturally tangent vector fields, e.g. u : 
Z) B ^ TB. In order to make use of exterior calculus, a corresponding displace- 
ment covector field u : E^ — > T* E^ is introduced, which is derived trivially from 
the displacement vector field by index lowering: 



Here, g{u) is the action of the covariant metric tensor on the vector u. The 
superscript b indicates the index lowering process or association of a covector field 
to a vector field, respectively. The inverse action is the index raising operation 



(2.1) 



/iAit + (A + Ai) V(V • -u) = 



(2.2) 



u :^ if' = g{u) 



u = u 
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2.2. Vector Differential Operators in Exterior Calculus. 

Lemma. For the differential operators grad, div and rot, applied on vector or scalar 
fields, respectively, the following relationships are valid: 



(2.3) V/ = df* 

(2.4) y -li = *d*u 

(2.5) V X IZ" = {*du)* 

Proofs for these relations can be found for example in |Fla90| . 

2.3. Cauchy Navier Equation. 

Theorem 1. Given a displacement vector field u that satisfies the Cauchy Navier 
eguations, then the corresponding covector field u satisfies 

(2.6) (A + 2^i)d5u + fiSdu = 

Proof. The correspondence between the divergence or rotation of vector fields and 
exterior derivatives of differential forms is well known. Since, however, there is no 
equivalent of the gradient of a vector field in exterior calculus, the Laplace operator 
has to be restated using the well known formula 

(2.7) Au = V • Vu = V(V • u) - V X V X u 

With this, the Cauchy Navier equation appears in the equivalent form 

(2.8) (A + 2/i)V(V • m) - X V X u = 

Now, using [^31 to 12. 51 we can restate in terms of the corresponding displace- 
ment one- form u = v}': 

(2.9) (A + 2^)d* * u - ^ * = 



A somewhat more concise notation of l2.9l uses the codifferential operator, defined 
in as 

(2.10) 5:^r{E^) n^'-^E^) 

-*d*uj ifujen^{E^)orujen^{E^) 

^ *d*u tfujen^{E^) 

With this we arrive at the proposed formula 



(2.11) 



(A + 2fi)dSu + fiSdu = 



□ 
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2.4. Boundary Conditions. The treatment of displacement boundary conditions 
are trivial: the displacement vector field on the boundary just has to be mapped 
to its corresponding covector field. 

Corollary. Given a displacement vector field u that satisfies the Cauchy Navier 
equations and assumes given values u\dB on the boundary dB, then the correspond- 
ing covector field u satisfuina \2.6\ assumes the values g{u) on the same boundary. 

Traction boundary conditions, however, are much harder to handle. It would 
be desirable to find a way to express traction forces or the elastic stress tensor in 
terms of the displacement covector field, resulting in boundary conditions that can 
be stated as a set of differential equations in the discplacement covector field. This 
will be the content of the next and largest part of this paper. 

3. Traction Boundary Conditions 

3.1. The Linear Strain Tensor. According to |Mar94| . the strain tensor G 
T2 (i?'^) of linear elasticity ^ = eij-e* (i) = ^(Vm + Vu"^) can also be written as 
the Lie derivative of the metric tensor with respect to the displacement vector field 
u: 



(3.1) ^ - ^inV 

with the metric tensor — gijdx^^dx^ . Beware that we need the discplacement 
vector field here, not the corresponding one-form. This is unfortunate, because we 
now need to map one-form into vector field components and vice versa. However, 
in the case of orthonormal bases and physical vector components this is simple: the 
components are unchanged. Otherwise, index lowering or raising has to be used. 

As an example, using a general canonical vector base, it shall be shown, that the 
formula above is equivalent to results of classical tensor analysis. 



Corollary. The components of \3.1\ satisfy 



(3.2) ^ = eijdx'' ® dx^ = ^{u'^gij.k + gkju'^i + gikUj)dx'' (g) dx^ 

Proof. This can be proven by direct calculation. 

The Lie derivative adheres to the product rule for tensor products: 

(3.3) Ls{T ® S") = LsT (g,S + T(g) L^S 

This holds for tensor fields of any rank, and therefore especially for scalar func- 
tions and one- forms. Thus we have 



(3.4) Laigijdx^ ® dx^) = {Lfigij)dx^ ® dx-' + gijL^dx^ (E) dx^ + gijdx'' ® L^dx^ 

Cartan's „magic formula" can be applied to calculate the lie derivatives of one- 
forms: 



(3.5) Lga — X ■ da + d{x ■ a) 

Applied to the canonical basis covectors dx^and afterwards to 13.41 
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(3.6) Lsdx^ = du' = u\kdx'' 

La{gijdx^ ® dx^) = u'^gij.kdx^ (g) dx^ + giju'j.dx^ ® dx^ + giju\dx^ (g) dx'^ 

(3.7) = iu''gij^k+gkju'^i + giku'^j)dx^ (S)dx^ 
Thus we have 

^ = eijdx^ (g) dx^ = ^(u'°5ij,fc + Qkju'^i + gikUj)dx'' (g) dx^ 

□ 

Corollary. The components of the linear elastic strain tensor defined as lij = 
\{ui\j + Uj\i) satisfy\3lE 

Proof. Let's first rewrite = ^{u^j + Uj^^) involving the ChristofFel symbols of 
the second kind: 

1 

2 

In the last step we used the symmetry property of the Christoffel symbol T^^ 

ji' 

Now we raise the index of the UiComponents and substitute T'^j — \g^™{gim,j 

= ^((Sifew'') J + {g'jku^),t) - ^g^'^igtm,] + .g^™,* - gtj,m)uk 

= ^{gik,ju'' + gikUj + g]k,iu'' + QjkU^i - gim.jU"^ - gjm.,iU"^ + 5ij,mW™) 

= \{gik,ju^ + ffifeWj + g]k,iu'' + gjku^i - gik,ju'' - gjk,iu'' + gij,ku^) 



(3.8) = -(liij + WjvO - Tj^j-Wfe 



J 

1 

A comparison shows that — e^j. □ 



(3-9) = oigikUj + gjku'^, + gij.ku'') 



3.2. The Stress Tensor. For isotropic homogeneous bodies, the stress and strain 
tensor of linear elasticity are related as follows |Mar94| : 



(3.10) V = AeV + 2Ai^ 

Here, e = V • u is the volume expansion under deformation. For a displacement 



one-form u we have e = —Su = *d * u. Since ^ = \L^^ we get: 



(3.11) W = -X6uY + l^LsY 

Beware that as before, unfortunately, we have to deal with the displacement 
one-form and vector field at the same time. 
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3.3. Traction Boundary Conditions. 

Lemma. Given a traction vector field t : dB 'D dS — >■ TB on apart of Bs boundary. 
The corresponding traction covector field t :=t^ then satisfies 



(3.12) t := P = ^ • n = Aen + • n 

Proof. This is a trivial implication of cauchy's stress theorem □ 

Theorem. Given a traction covector field t : dB D dS —^T*B and a displacement 
covector u field satisfvina \2.b\ The thus given traction boundary conditions are sat- 
isfied, if the displacement covector field and its corresponding vector field u adheres 
also to the following equation on dS: 

(3.13) t = —XSun + ii{d{u ■ n) + u ■ dn + [ft, u]'') 

where n is a normal vector at a given point and n — ffi its corresponding one-form. 

Proof. Since the Lie derivative is also a derivative with respect to scalar products 
{Ls{T ■ S) — LgT ■ S + T ■ LgS for any two tensors T and S) we can write: 

2^ • n — Lzlf ■ n — Lz{^ ■ n) — ^ ■ L^n 

= L^n — [u, fi]^ — d{u • n) + u • (in + [ft, u]'' 

Now only the part Xen is missing from 13.121 which is equal to —Xdun. □ 

Now let's assume we have a coordinate system on which one of the canonical base 
vectors equals the surface normal ft. In accordance with cylindrical or spherical 
coordinates, let's call this coordinate function r, the canonical base vector along 
this coordinate dr and the canonical base covector dr. Let's also assume that n 
(and thus n) are normalized, as well as dr and dr. In this case ft — dr and n = dr. 
Then we can reformulate 13.131 according to the following 

Corollary. Given a coordinate system with a coordinate function r satisfying the 
description above, a traction covector field t : dB D dS — ?• T*B and a displacement 
covector u field satisfying \2.6\ the latter have to adhere to the following equation 



(3.14) t = -{XSu)dr + fi{du'' + [dr,u]'') 

Proof. We look at all the substitions one after the other. ~X5un = —{X5u)dr just 
by the definition n = dr given above. 

d{il ■ n) = du^ follows from the fact that u ■ n = u ■ dr ~ . 

u ■ dn vanishes because dn — ddr = 0. 

[ft, u]^ — [dr, u]^ also by definition. □ 
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4. Conclusion 

In the first part of this article we have derived Cauchy Navicr's equation of static 
equilibrium in linear elasticity using only differential forms instead of vector fields. 
The second part treats the harder problem of expressing strain and stress tensors 
in a more modern differential geometric way, and, finally, how to deal with traction 
boundary conditions. 
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